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(54) Interpolation filter 

(57) An interpolator for a discrete time signal is pro- 
vided, which performs the steps of identifying a plurality 
of interpolation points which interpolation points are 
symmetrically temporarily displaced about a central ref- 
erence point of a plurality of discrete time samples of a 
signal for which interpolated signal samples are to be 
calculated at the interpolation points, determining a plu- 
rality of impulse response coefficients appertaining to 
the interpolating function at each corresponding sam- 
pling time of the said plurality of signal samples, and 
contemporaneously calculating a plurality of interpo- 
lated discrete time signal samples corresponding to the 
interpolation points in combination with a pre-addition of 
at least one pair of discrete time signal samples. 
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mated a, a particular point from vSu s ^SS^^S^J^ V" 6 * 3 *" clon ° r si 9 nal is ^ 
ing. interpolation is usedto increase a sample S SSKSSL?^ . 6 COnteXt 01 di9ita ' Si9nal process - 
time signal samples representative of the tiS5a£ £££ J, ? V 11,716 Si9 " al Samp,es " 77,6 discr ete 

sampling interval of thesignal. The uniSZ^^SS^S^^ h accordan « with a 

Placements. Interpolated signaisampTe are ^aLTc^„?f ^ S ' 9na ' Samp ' eS at '"^mediate time dis- 

re9 T;Tn^ 

"^Z^ or indeed of any Pinite 

A characteristic of a linear phase impulse response ZT^t 1 8 l,near phase is known as 'Pre-addition'. 
tribution about a centre or^ere^ntT^a^^T.I "T^** reSP ° nSe have a ^etrica. dis- 
response which is characterised in that each Hr^u.s^nSllr T t ^ ttMry °* the ' inear phase im P ulsa 
in an equiva.ent temporal position aboul a 3fe 0 ^ e f en ^ 

response a , = a ., The pre-addition structure ierati to ^L?!!^ ^ ^ ** 8 ' inear phase im P ulse 
of discrete time signal samples held within th?S TI^/pbm ^T"^ inter ?° ,ator by summing pairs 
response coefficients wherein the pairs of signal samS^rr^L 'nterpolator. before multiplication by impulse 
cients of the interpolator, or. in other words ^'^^0^ Symmetrical ^ <* 'mpu.se response com- 

struS^gl™ interPO^ng fi,te, with a pre-addition 

t-me signal samples via a conductor 9. Also^own in FiouTe it ?J5 2 r 9 ^ 1 ' ?" 3 ' 4> 5> 6 " 7 ' 8 " ,ed with discrete 
which is representative of a discrete reference cTeTnnS f f °u " n6S 1 °' 1 1 ' 12 ' 13 ' 14 ' 15 - 16 - 17 - each of 
3. 4. 5. 6. 7. 8. of the shift register r taXf SEETtf TS S Tf ^ ° f "* ° f the StageS 1 ' 2 ' 

represented as Xj where i represents the relafivePoSn in i ottho " ' 6 ' 1? ' ,s sh °wn with a signal sample 
centra, or reference point x 0 . Also shown in RgureTt a co^ce^alTneT ^T* ° f the With respect to a 
The interpolating function 18. is shown to be Symmetrel JSilf " rep | esentat,ve of an interpolating function, 
to a temporal position of a interpolated signafSe to b e calLlt^T P ° Si,i ° n y °' whicn <*>"esponds 

function 1 8. each of the coefficients rTpresenS AS 3 feSU,t ° f the symmetry of interpolating 

as aforementioned in that a 1 = * a^T^ s £Z£ ST* ° f tha interpo,a ^ Action is symmetrical 
Plurality of summing amplifiers 19 20. 21 22 whS, SeraS totl l i = ^ • Connected «° *e shift register V is a 
corresponding coefficients of the impulse r'esoons TcZeM^JZ^r^" ° f diSCrete fime Si 9 nal sam P les «* 
At an output of each of the summing ampS 19 20 2,2 ^T^T 18 whichare ^metrical in that a i= a, 
sponding composite coefficient a ( of the inte Satina tl hT,? ? $Um * subse ^"ently scaled by corre- 
Each output from the multipliers 23 24 25 £° X ™ ?" ? * 3 Corres P° n ding plurality of multipliers 23. 24. 25. 26. 

•^?^^^o^^2* ^"-pSSSL, filter shown in Figure 1 
achieved by pre-adding J^^^^^^J^ST^ ^ *** ** interpo,ati ° n - ™s Jving is 
oal impu.se response coefficients a^d mu^ 

impulse response coefficients. The result crfTSiSS«i^ ^ k^ 1 " by ° nG ° f the pair of ^spending 
after summed to form an interpolated signal s^pTe P "* P3 ' r 0< diSCre,e time si 9 nal sam P'es are there 

desc^^ with a pre-addition structure as hereinbefore 

structure is conditioned on the irrpuTse re S pote of SESSJST? ^ Pr ° ViSi ° n ° f a pre - additio " 

possessing a necessary symmetry. Although an iS5 , N 3 BnMr PhaSS reSponse and hereby 
Phase, an impulse response representative o" a JSl^S^ ^«**« filtar exhibit a lineaV 
about an interpolation point of the funSon coneSSc ^to aTel^,' nter P 0lat,on n«y not be symmetrical 

the use of a pre-addition structure co n-espond.ng to a temporal posrt,on of a interpolated sample, preventing 
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steps of; 

(i) identifying a plurality of interpolation points which interpolation points are substantially equally temporarily dis- 
placed with respect to at least one point of a plurality of discrete time samples of a signal for which interpolated sig- 

s nal samples are to be calculated at the said interpolation points; 

(ii) for each of the interpolation points, determining a plurality of impulse response coefficients appertaining to the 
interpolating function at each corresponding sampling time of the said plurality of signal samples, and 

(iii) contemporaneously calculating the interpolated signal samples corresponding to the interpolation points in 
combination with sum and difference operations of at least one pair of impulse response coefficients, and at least 

io one pair of discrete time signal samples. 

By selecting a plurality of interpolating points, each of which is displaced in time by an equal amount from a central 
reference point, the calculation of interpolated signal samples corresponding to the interpolation points may be sub- 
stantially reduced by utilising a relative symmetry of the impulse response coefficients associated with respective inter- 
15 polation points. 

One embodiment of the present invention will now be described by way of example only, with reference to the 
accompanying drawings wherein; 

FIGURE 2 is a representation of two versions of an interpolation function with associated coefficients correspond- 
20 ing to two interpolating points for an odd order interpolator. 

FIGURE 3 is a representation of two versions of an interpolation function with associated coefficients correspond- 
ing to two interpolating points, for an even order interpolator, and 
FIGURE 4 is a schematic block diagram of an even order interpolating filter. 

25 A Finite Impulse Response Filter (FIR) which operates to convolve an impulse response of the filter with a discrete 
time signal may be an even order filter or an odd order filter. An even order filter comprises an impulse response with 
an even number of impulse response coefficients whereas an odd order filter embodies an impulse response with an 
odd number of irrpulse coefficients. An odd order impulse response is characterised by a fact that a discrete time sam- 
ple will tie at a central point of the impulse response of the filter. An even order impulse response is characterised in that 

30 a sample does not lie at a central point of the impulse response of the filter. 

An illustrative representation of an interpolation function in continuous time form is shown in Figure 2. In Figure 2 
two versions of the interpolating function are shown 29a, 29b, which correspond to a calculation of two interpolated sig- 
nal samples y 0 y v Also shown in Figure 2, are a plurality of discrete time sampling points represented by lines which 
are labelled with x_j where index i is in the range -3 to 3. In this case an interpolator for interpolating the signal sample 

35 is an odd order interpolator. At each of the discrete time sampling points x_ 3 , x_ 2 . x^ , xq, x t , x 2 , x 3 , a corresponding sam- 
ple of the interpolating functions 29a, 29b is shown by impulse response coefficient a ; where i is also in the range -3 to 
3. The impulse response coefficients, a { form an impulse response corresponding to a discrete time representation of 
the interpolating function 29a, 29b which serve to form interpolated samples y 0 , y 1 . As a result of the fact that the inter- 
polated sampling points Vj y 1 are temporarily displaced with respect to the discrete time sample instant xo by an equal 

40 amount, the impulse response coefficient for each of the interpolating points ,y 0 are substantially the same, in that 
the coefficient aj where i is -3 to 3 are the same for a calculation of an interpolating sample at the interpolation point y 0 
as they are for the interpolating point Y v The coefficients for calculation of respective interpolated samples y 0 yi are 
therefore symmetrical. It is this symmetry between calculation of interpolated samples that provides a saving in compu- 
tation of the interpolated samples at points y 0 y 1( in accordance with the following method. 

45 In the example shown in Figure 2, the interpolation point y 0 corresponds to a temporal position, which is a quarter 

of the sampling interval earlier than a current input signal sample position corresponding to x 0 . 

The interpolator point y-i corresponds to a temporal position, which is a quarter of the sampling interval later than 
the current input signal sample corresponding to x^ The coefficients for the interpolation point are shown on the dia- 
gram as aj with/c{-3...3). Thus 

50 

3 
;»-3 

55 

and, more generally for an odd order interpolator :- 
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Xo- Z a < x i 
/—A 



5 S^£^~~is^'S£==»''—- 



/=-3 



/5 



and, more generally:- 



/=-* 



We may express the above as follows:- 
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>o = *o*o + 2, a iXi + X = OoX 0 + + ^ J 

and 

i k k 

-v, = *o*o + X fl _ iJf , + = a 0 x 0 + + <V«L,) 



30 Now form 



(1) 



1=1 



where a'- 



/ = 0 — 

2 



and 



t >0 a±£=i 



(2) 

ft - -Vn-y. _ 1 V/ 

£ = £2>, - "-,)(*, - AT.,) = £^ ~ *-.) 
'=1 ._■ 
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where q."= °' °- 
2 



(1) 



(2) 



The interpolated signal samples corresponding to the sampling 



points y 0 . y, are thereafter calculated in accord- 
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ance with the sum of a and p for y 0 , and a difference between a and p for as represented by equations (5) and (6). 

y 0 = a + p (5) 

5 y^a-p (6) 

A representation of two versions of the interpolation function for calculation of two interpolated samples at positions 
y 0 y-| for an even order impulse response is shown in Figure 3 where parts or elements also represented in Figure 2 
bear identical numerical or alpha numeric designations. In Figure 3, the interpolating function 29a, 29b, are shown to 
w be sampled in accordance with a plurality of discrete time sampling positions xi as for Figure 2, but unlike Figure 2, there 
are now eight sampling positions and correspondingly the impulse coefficients a, for each of these even number of sam- 
pling positions are adjusted in accordance with the temporal position of each of the discrete time sampling points xi 
where now i is in the range -4 to +3. In this case we have 

15 K-1 

y 0 = Z a i x i 

/=-/< 

20 for a filter of order 2k{k - 3 in Figure 3). With a symmetry similar to that described earlier, we have 

= Z a -iM)*r 

25 

We can express the above as> 



-1 K-1 /C-1 

Vo = Z a / x / + Z a <*/= Z < a / x < + a -</ + i) x -</+i)) 

i=-k /«0 /=0 



30 



and 

35 

-1 /c-1 *-1 

y-i - Z a -(/+i) x / + Z a -(/+i) x / = Z ( a -(/+i) x / + a / x -(/+i>) 

/—ft foO /=0 

40 

Now form 

a = ° 2 1 = 2 Z ( a / x / + a -(/+i) x -(/+i) + a «(/+i) x / + a / x -(/+l)) ( 3 ) 

ft 



<* = g Z < a / + a -(/ + i)H x / + *-(/+d> = Z a /( x / + x -(/ + D ) 



/-1 /01 



50 

where 



55 a / = 2^ 



and 
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yo-y* 1 H 



U+i)) (4) 



1 * * 

p - I Z < a < • a -(/ + i))(^/ " * ( , + d = 2 a 0 " (x, - x. ( ) 
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where 



" a i' a -(/+1) 
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15 muXttnfea^ ^ ~ ^ — ■> - * — * where a and p are generated with only * 

« we^rform WhiCh US6S th6 ** 6XpreSsions to aerate a and p. For 

the taps a,:, / e{0 k^^^^^^^TT^^ pre " addin 9 the centre to itself) and multiply by 

pies (in this case the ^ZZ^*™ «" F °' P we perfo ™ a P^subtractJon over the sam- 

^.MI.^.ThuswepedoTm^ 

p. Given a and p we can easily derive vn and v <;inr 0 „ . ^ P multiplications to generate a and 

have generated both y 0 and yXZ^ooTL^ mu.«l»« ^ = a - p as hereinbefore explained. Thus we 

33 £ 7s U t6 4 37 d ir^? 9na] t ,6d ,r ° m 3 C ° ndUCt0r 30 ' t0 8 shift re 9 ister ». comprised of eight stages 31 32 
tractor3 3 9 5 'i conned t^^ 
connactedtoasecondst^eSandto^ 

41. is connected to the third staoe 33 an^ Jl °«J ' L . u re9 ' ster58 ' ""areas ttie third adder subtracter 
tractor 42. is connected to the fouSa Q ?34 and o thS ♦ ' "2? ? ? r " BW-r ^ ^ the fourth adder su °" 
four adder subtracters 39. 40 4^ T2l s ?es P ectiv2 35 ' ° f f3 Shm regiSter 58 An output of each °' tha 
An output of each of the four mulipLs 2 46 ,s ^SSSS . ^ ^ ° f *"* ^ 45 and 46 
47. is connected to one positive in^^^^ ° f an adder 47 An «*« * *• adder 
negative input of a second arithmetic uniM9 The oT« f IS S T ' 6 4?> ' S alS ° COnnected to a 
An output of the latch 50 is connected ToTJ™^ ♦ ? a ' S ° connected t0 an input of a latch 50. 
input oTthe second arithmetic unZt A ctocK inou^heS^ " ' *™ C 48 and 3 second 

is also connected to a clock ir^ «f££ oT^o D tvl P »J T"£Tr? 3 m " tBr C,0Ck Via 3 COnduCtor 51 ■ which 
to a clock input of the shift renter 58 and reso-ti^, T P ^ A ° ° UtpUt ° f the D - ,ype ,,ip - ,lo P is connected 

42. The Q output from the D 2 ,f c ill , 2? Y * °J C ° ntr °' ' npUt ° f each of 106 adder subtracters 39. 40. 41. 

An output fronTeach of the ses 52 S V^T^^ ^ **« eaCh ° f ,0ur switches «. 53. 54 55. 

43. 44. 45. 46. ' ' ^ " reS P eCt,vel V reeled to a second input Qf egch Qf ^ mu|tjp|jefS 

In Figure 4. the shift register 58. clock in signals samples x where si Th= ■ _. 

output for / = 0 and the index increases fram rSJSTiT^ 1 > ™ S ri9htm ° St adder ^btractor provides the 

adder/subtracter is J^^SSJSS SlZlEH; ? 0U,PUt iS mU,tip,ied * either a i < if *• 
of the switches. ' " me aader/subtractor 15 subtractmg) ■ The cho,ce of a, or a ; is according to the state 

During the interval between successive cvcles of thp ctamnio 
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25 



a The second phase begins with the next rising edge of the master clock. This clocks the value of a into the latch ^5 7 
and sets the sample clock low. In this state the circuit is arranged to compute the value of Once the «^d«be 
of p becomes valid (towards the end of the master clock cycle), the outputs of the adder and ™™*<*°<V™** 
y 0 = a + p and y, = a - p respectively. Thus, in this implementation, the two outputs are generated contemporane- 

° USl As will be appreciated by those skilled in the art various modifications may be made to the arrangement of the inter- 
polating filter without departing from the scope of the invention, in particular it may be appropriate to interpolate severe J 
values between the regular input samples. There may. nevertheless, remain a symmetry between the coefficients 
needed to generate half of a subset of these samples and the other half of that subset. The pr.nc.ples hereinbefore 
,o described may therefore be extended to this case by generating the same filtering operations to generate each suitable 
pair of interpolated values. Note, however, that a further saving is possible in this case because the required Po- 
tion and pre-subtraction operations are common to all required interpolated values. Thus these operates need only 
be performed once. 

is Claims 

1 . A method of interpolating a discrete time signal comprising steps of; 

(i) identifying a plurality of interpolation points which interpolation points are substantially equally temporarily 
displaced with respect to at least one of a plurality of discrete time samples of a signal for which .nterpolated 
signal samples are to be calculated at the said interpolation points; 

(ii) for each of the interpolation points, determining a plurality of impulse response coefficients appertaining to 
the interpolating function at each corresponding sampling time of the said plurality of s.gnal samples, and 
(Hi) contemporaneously calculating the interpolated signal samples corresponding to the .nterpolat.on po.nts in 
combination with sum and difference operations of at least one pair of impulse response coefficients, and at 
least one pair of discrete time signal samples. 

2. A method of interpolating a discrete time signal as claimed in Claim 1. wherein the said plurality of impulse 
response coefficients of the interpolation function determined in step (ii) of the said method are substantially the 

30 same for each of the said interpolation points. 

3. A method of interpolating a discrete time signal as claimed in Claim 1 and 2. wherein step (iii) of the said method 
comprises the steps of; 

(iv) forming first intermediate data a, in accordance with at least one summed pair of discrete time signal sam- 
ples in combination with at least one summed pair of impulse response coefficients; 

(v) forming second intermediate data p.in accordance with a difference between at least one par of discrete 
time signal samples, in combination with a difference between at least one pair of impulse response coeffi- 
cicnts snd 

(vi) calculating the interpolated signal samples in accordance with a difference between the first intermediate 
data a, and the second intermediate data p. and a sum of the f irst intermediate data a. and the second inter- 
mediate data p. 

A method of interpolating a discrete time signal as claimed in Claim 3, wherein each signal sample of the said pair 
of at least one pair of signal samples in step (iv) of the said method for forming the first intermed.ate sum a, is sub- 
stantially equally temporally displaced with respect to the central reference point, and the said at least one pair of 
impulse response coefficients are representative of samples of the interpolation function at temporal positions 
appertaining the temporal displacement of the said at least one pair of signal samples with respect to the central 
reference point. 

A method of interpolating a discrete time signal as claimed in Claim 3 or 4, wherein each signal sample of the said 
pair of at least one pair of signal samples in step (v) of the said method for forming the second intermediate sum p. 
is substantially equally temporally displaced with respect to the central reference point, and the sad at least one 
pair of impulse response coefficients are representative of samples of the interpolation function at temporal posi- 
tions appertaining the temporal displacement of the said at least one pair of signal samples with respect to the cen- 
tral reference point. 

A method of interpolating a discrete time signal as claimed in Claim 3 to 5. wherein each of the said pairs of signal 
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7. 



8. 



same as *, saidpairs „ rCtlZZSSZ T"?" SUm "• are "-"SfSl 

mediate sum p. ^ »» rr°e"iO'ems in step (v) of the said method tor forming the second inter- 

^stro^s 

1 * * 
a - * o*o + 2 -I (a, + a .,)(*,. + * .) = J */'(*, + x , 

/«<> 



wherein 







f'- 


2 


\i>0 






2 



accordance with equation (2) as follows; 9 the SeCOnd ,nterm ediate sum 3, is calculated in 

where 

« - 1 2 (a, . ».,,.„)„, «. , w , w 

Is1 * » 



where 



rsi^aTJa 3 ^ the interpolation function, and ^represents the i- 

accordance with equation (4] las follows 9 S6COnd interme diate sum 0, is ca.cu.atea i 
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s where 

_" a /" a -(/fi) 
a i = 2 

10. A method of interpolating a discrete time signal as claimed in Claims 8 or 9, wherein the plurality of interpolation 
w points are two interpolation points, designated y 0 and y 1f and wherein y 0 = a + p and y t = a - p . 

1 1 . An interpolating filter for generating a plurality of interpolated signal samples of a discrete time signal, comprises a 
" shift register, which shift register serves to store samples of the discrete time signal fed thereto, a plurality of switch- 
able adder/subtractors connected to the said shift register, each of which said adder/subtractors operates in a first 

is mode to add a pair of signal samples communicated from the shift register, and in a second mode to subtract the 
pair of signal samples from each other, and which adder/subtractors are switched between first and second modes 
in accordance with a control signal, a plurality of scaling means, each of which scaling means operates to scale 
signals communicated from an output of the said adder/subtractors with data communicated from one of a plurality 
of controllable switches, the data fed from each controllable switch being either a sum of two interpolating coeff i- 

20 cients in a first position of the said switch, or a difference between the interpolating coefficients in a second position 
of the said switch, which position of the switch is controlled in accordance with the said control signal, a summer 
connected to an output of each of the plurality of scaling means, which summer operates to sum scaled outputs 
from the said adder/subtractors, and an interpolating controller which operates to generate the control signal and 
the interpolated signal samples in accordance with sum and difference operations on signals provided at an output 

25 of the summer. 

12. An interpolating filter as hereinbefore described with reference to Figures 2, 3, and 4. 
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